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Abstract 

The modified theories of gravity, especiahy the f{R) gravity, have 
attracted much attention in the last decade. In this context, we study 
the exact vacuum solutions of Bianchi type /, /// and Kantowski- 
Sachs spacetimes in the metric version of f{R) gravity. The field 
equations are solved by taking expansion scalar 9 proportional to shear 
scalar a which gives A = B"^, where A and B are the metric coeffi- 
cients. The physical behavior of the solutions has been discussed using 
some physical quantities. Also, the function of the Ricci scalar is eval- 
uated in each case. 

Keywords: Bianchi type /, Bianchi type ///, Kantowski-Sachs and f{R) 
gravity. 
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1 Introduction 

The late time accelerated expansion of the universe has attracted much at- 
tention in the recent years. Direct evidence of cosmic acceleration comes from 
high red-shift supernova experiments [1]. Some other observations, such as 
cosmic microwave background fluctuations p] and large scale structure [3], 
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provide an indirect evidence. These observations seem to change the entire 
picture of our matter filled universe. It is now believed that most part of the 
universe contains dark matter and dark energy. The modifications of general 
relativity (GR) seem attractive to explain late time acceleration and dark 
energy. 

Among the various modifications of GR, the f{R) theory of gravity is 
treated most seriously during the last decade. It provides a natural gravita- 
tional alternative to dark energy. It has been suggested that cosmic acceler- 
ation can be achieved by replacing the Einstein-Hilbert action of GR with a 
general function Ricci scalar f{R)- f{R) theory of gravity has been shown [1] 
equivalent to scalar-tensor theory of gravity that is incompatible with solar 
system tests of GR, as long as the scalar field propagates over solar system 
scales. The explanation of cosmic acceleration is obtained just by introduc- 
ing the term 1/R which is essential at small curvatures. Capozziello et al. 
[5] have shown that dust matter and dark energy phases can be achieved 
by the exact solution derived from a power law f{R) cosmological model. 
The f{R) theory of gravity is considered most suitable due to cosmologically 
important f{R) models. These models consist of higher order curvature in- 
variants as functions of the Ricci scalar. Viable f{R) gravity models ^ have 
been proposed which show the unification of early-time inflation and late- 
time acceleration. The problem of dark matter can also be addressed by 
using viable f{R) gravity models. There are some other useful aspects [7] of 
f{R) gravity. It gives an easy unification of early time inflation and late time 
acceleration. It can be used for the explanation of hierarchy problem in high 
energy physics. It also describes the transition phase of the universe from 
deceleration to acceleration. Thus f{R) theory of gravity seems attractive 
and a reasonable amount of work has been done in different contexts. 

Lobo and Oliveira |8j constructed wormhole geometries in the context 
of f{R) theories of gravity. Cognola et al. [9] investigated f{R) gravity at 
one-loop level in de-Sitter universe. It was found that one-loop effective ac- 
tion can be useful for the study of constant curvature black hole nucleation 
rate. Multamaki and Vilja [TD] investigated spherically symmetric vacuum 
solutions in f{R) theory. The same authors [H] also studied the perfect 
fluid solutions and showed that pressure and density did not uniquely deter- 
mine f{R). In a recent paper, Sharif and Kausar |T2] studied no n- vacuum 
static spherically symmetric solutions in this theory. Capozziello et al. [13] 
explored spherically symmetric solutions of f{R) theories of gravity via the 
Noether symmetry approach. HoUenstein and Lobo analyzed exact so- 
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lutions of static spherically symmetric spacetimes in f{R) gravity coupled to 
non-linear electrodynamics. Azadi et al. [TB] studied cylindrically symmetric 
vacuum solutions in this theory. Momeni and Gholizade [16] extended cylin- 
drically symmetric solutions in a more general way. Reboucas and Santos [17] 
studied Godel-type universes in f(R) gravity. We have explored static plane 
symmetric vacuum solutions [IB] in f{R) gravity. The field equations are 
solved using the assumption of constant scalar curvature which may be zero 
or non-zero. In a recent paper, Babichev and Langlois fT9] studied relativistic 
stars in this theory. 

Friedmann- Robert son- Walker (FRW) models, being spatially homoge- 
neous and isotropic in nature, are best for the representation of the large 
scale structure of the present universe. However, it is believed that the 
early universe may not have been exactly uniform. Thus, the models with 
anisotropic background are the most suitable to describe the early stages 
of the universe. Bianchi type models are among the simplest models with 
anisotropic background. Many authors [20]- [27] explored Bianchi type space- 
times in different contexts. Kumar and Singh [2H] studied solutions of the 
field equations in the presence of perfect fiuid using Bianchi type / space- 
time in GR. Moussiaux et al. [29] investigated the exact solution for vacuum 
Bianchi type-/// model with a cosmological constant. Lorenz-Petzold [20] 
studied exact Bianchi type-/// solutions in the presence of electromagnetic 
field. Xing-Xiang [30] discussed Bianchi type /// string cosmology with 
bulk viscosity. He assumed that the expansion scalar is proportional to the 
shear scalar to derive the solutions. Wang [31j investigated string cosmo- 
logical models with bulk viscosity in Kantowski-Sachs spacetime. Upadhaya 
[32] explored some magnetized Bianchi type-/// massive string cosmological 
models in GR. Recently, Hellaby [53] presented an overview of some recent 
developments in inhomogeneous models and it was concluded that the uni- 
verse is inhomogeneous on many scales. 

The investigation of Bianchi type models in alternative or modified theo- 
ries of gravity is also an interesting discussion. Kumar and Singh [31] inves- 
tigated perfect fiuid solutions using Bianchi type / spacetime in scalar-tensor 
theory. Singh et al. [3S] studied some Bianchi type-/// cosmological mod- 
els in scalar-tensor theory. Adhav et al. [36] obtained an exact solution of 
the vacuum Brans-Dicke field equations for the metric tensor of a spatially 
homogeneous and anisotropic model. Paul et al. [37] investigated FRW cos- 
mologies in /(/?) gravity. Recently, we [33 139] have studied the solutions of 
Bianchi types / and V spacetimes in the framework of f{R) gravity. 
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In this paper, we focuss our attention to explore the vacuum solutions of 
Bianchi types /, /// and Kantowski- Sachs spacetimes in metric f{R) gravity. 
The field equations are solved by taking expansion scalar 6 proportional to 
shear scalar a which gives A = 5", where A and B are the metric coefficients. 
The paper is organized as follows: A brief introduction of the field equations 
in metric version of f{R) gravity is given in section 2. In section 3, 4 and 5, 
the solutions of the field equations for Bianchi types J, III and Kantowski- 
Sachs spacetimes are found. Some physical parameters and the functions of 
Ricci scalar are also evaluated in the context of these solutions. In the last 
section, we discuss the results. 

2 Some Basics of f{R) Gravity 

The metric tensor plays an important role in GR. The dependence of Levi- 
Civita connection on the metric tensor is one of the main properties of GR. 
However, if we allow torsion in the theory, then the connection no longer 
remains the Levi-Civita connection and the dependence of connection on the 
metric tensor vanishes. This is the main idea behind different approaches of 
f{R) theories of gravity. 

When the connection is the Levi-Civita connection, we get metric f{R) 
gravity. In this approach, we take variation of the action with respect to the 
metric tensor only. The action for f{R) gravity is given by 

j V^{f{R) + Ljd\ (1) 

where f{R) is a general function of the Ricci scalar and is the matter 
Lagrangian. The field equations resulting from this action are 

F{R)R^, - \f{R)g,. - V^V,F{R) + g,.nF{R) = T^, (2) 

where F{R) = df{R)/dR, □ = V^V/^, is the covariant derivative and 
is the standard minimally coupled stress energy tensor derived from the 
Lagrangian L^. Now contracting the field equations, it follows that 

F{R)R - 2f{R) + 3nF{R) = T 

In vacuum, this reduces to 

F{R)R - 2f{R) + 3nF{R) - 0. (3) 
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which imphes that 

This gives an important relationship between f{R) and F{R) which will be 
used to simplify the field equations and to evaluate f{R). 
Now we consider the metric 



ds"" = de - A\t)dr^ - B\t)[de^ + l\e)dct)\ (5) 



where 



9'^ when k=0 (Bianchi I model), 

t^{9) = ^ sinh^ when k=-l (Bianchi III model), 

sin^ ^ when k=l (Kantowski-Sachs model). 

Here k is the spatial curvature index and the above three models are 
Euclidian, semi-closed and closed respectively. 



3 Bianchi Type / Solution 

Here we shall find exact solutions of Bianchi / spacetime in f{R) gravity. For 
the sake of simplicity, we take the vacuum field equations. The line element 
of Bianchi type I spacetime is given by 

ds^ = dt" - A\t)dr^ - B\t)[d9^ + 9^d(P% (6) 

where A and B are cosmic scale factors. The corresponding Ricci scalar is 
given by 

^^A 2B 2Ab B\ _ 
^ = -2^1 + ^ + ^ + ^]' 
where dot denotes derivative with respect to t. Using Eq.(jl]), the vacuum 
field equations take the form. 



FiR)R^, - V^V,F{R) _ F{R)R - DF{R) 
9^v 4 



(8) 



One can view Eq.® as the set of differential equations for F{t), A and B. 
It follows from Eq.(|8]) that the combination 

^ F{R)R,,-V,V,F{R) ^ 
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is independent of the index and hence Afj_ — A^y = for all jj and u. Thus 
Aq — Ai = gives 

2B 2AB Af F 

Also, Ao- A2 = yields 

A B AB B^ BF F 

Now we give definition of some physical quantities before solving these equa- 
tions. 

The average scale factor a and the volume scale factor V are defined as 



a = VAB^, V = a^ = AB\ (12) 
The average Hubble parameter H is given in the form 

1 J 2B, , , 

The expansion scalar Q and shear scalar a are defined as follows 



" = "■'• = 1 + 4 



where ^ ^ 

^iiv = ^iUi,;aK + Uu;ah1l) " ^dh^^, (16) 



^Aii/ = 9tMu — Ufj,Uy is the projection tensor while = ^^(1, 0,0,0) is the 
four-velocity in co-moving coordinates. 

We have two differential equations given by Eqs. fllOflip with three un- 
knowns namely A, B and F. Thus we need one additional constraint to 
solve these equations. We use a physical condition that expansion scalar is 
proportional to shear scalar a which gives 

A = B"". (17) 
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Using this condition, Eqs. fllOflip take the form 

2B B"^ BF F 
-_ + 2„^ + „— -- = 0, (18) 

. B , r, ^ B BF F , . 

(^ + l);^ + (^^-2n-l)-^- — + - = 0. (19) 



Adding these, we get 



B , BF 

5 + (n + l)p + ^ = 0. (20) 



We solve this equation using power law relation between F and a [38] . 

F = kaJ^, (21) 
where k is the constant of proportionality, m is any integer and a is given by 

a = B— (22) 

Thus for m = 3, we obtain 

F = kB''^'^. (23) 
Using this in Eq. fl20l) . it follows that 

| + (2n + 3)|^ = 0. (24) 

Put B = g{B) in this equation, we get 

dg'^ ^ 4n + 6 2 



which leads to the solution 

9^ = -^^^ (26) 
where Ci is an integration constant. Hence the solution becomes 

ds' = A'dB' - B'^dr^ - B\de' + e^d<P^), (27) 
dB 

which can be written as 

ds^ = ^B^'^+^'dB^ - B^^'dr^ - B\de^ + e^d<p^), (28) 
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and after the transformations B = T, r = R, 6 = Q and = $, it takes the 
form 

ds^ = -T^^+^dT^ - T2"rfi?2 - T\d&^ + e^d^^). (29) 

The average Hubble parameter becomes 

_ + 2) 

while the volume scale factor turns out to be 

V = T"+2. (31) 

The expansion scalar 6 is given by 



_ y^(n + 2) 

while the shear scalar a becomes 

2 _ ci(n- 1)^ 

Moreover, the function of Ricci scalar, f{R), can be found by using Eq.(|l]) 

m) ^ (34) 

It follows from Eq.(I7j) that 

R = 2ci {n^ + 6n + 5)r-^"-l (35) 
Thus f{R) can be written as a function of R only 

/(i?) = -[2ci(n2 + 6n + 5)]^i?i (36) 

4 Bianchi Type /// Solution 

Here we shall find the solution of the Bianchi type III spacetime in f{R) 
gravity for the vacuum field equations. The line element of Bianchi type /// 
spacetime is given by 

ds^ = de - A\t)dr^ - B\t)[d9^ + sinh^ 9d(^% (37) 
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where A and B are cosmic scale factors. The Ricci scalar for this spacetime 
is given by 

A 2B 2Ab 1 B\ 

^a^^^^b'b^^W^ 

where dot denotes derivative with respect to t. Using Eq.([n]), the vacuum 
field equations take the form, 



R = -2[t + — + -rw-ir2+ml (38) 



2B 2AB AF F 

A B AB 1 B^ BF F _ 
~A~B^AB~B^^B^^W^~F~^' 

Here we also need one additional constraint to solve these equations. Using 
the same physical condition that expansion scalar 6 is proportional to shear 
scalar a, we obtain 

2B B^ BF F 
__ + 2„_ + „___ = 0, (41) 

(„ + l)| + („=-2„-l)|-§ + f ^ -i,. (42) 

Adding these, we get 

B , ,B^ BF 1 

B + (" + l)ij5 + ljF = -(;r3T)l55- («) 

Using power law relation between F and a, we obtain 

F = kB''+^. (44) 
Thus we get a differential equation with one unknown, 

| + (2„ + 3)|; = -^^^. (46) 



Put B = g{B) in this equation, we get 

dg'^ ^ 4n + 6 2 



dB + — « = - (^73T)B- («) 
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which leads to the solution 

,2 = ^ I (47) 

y 54n+6 (n-l)(2n + 3)' ^ ' 

where C2 is an integration constant. Hence the solution becomes 

ds^ = { _^ 1 )dB^ - B^^'dr^ - B^{de^ + sinh' ^0^), (48) 

_B4"+6 ~ (n-l)(2n+3) 

which takes the form 

ds^ = {— i ^ )dT^ - T^'^dR^ - T\dQ^ + sinh' e^^). (49) 

r4".+6 ~ (n-l){2n+3) 

where B = T, r = R, 6 = Q and = $. 

The average Hubble parameter becomes here 

n + 2 C2 1 ,1 ^ . 

3 ^i2^4n+8 (n- l)(2n + 3)T2J'' ^ 

while the volume scale factor turns out to be same as for the Bianchi type / 
spacetime. The expansion scalar 6 is given by 

'^-(" + 2)lf^- (,,_i)(i + 3)r^ l*. (81) 

while the shear scalar a becomes 

= - ')'[tS^ - (n-l)(2n + 3)T^ ^- ^''^ 
Moreover, the function of Ricci scalar, f{R), can be found by using Eq.(|l]) 

It follows from Eq.(^ that 

„ ^.C2(n2 + 6n + 5) 1 .Sn^ + 2n - 2, ^ 

which clearly indicates that f{R) cannot be explicitly written in terms of R. 
However, for a special case when n = —5, f{R) turns out to be 

^ ^ 12^3 
This gives f{R) only as a function of R. 
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5 Kant owski- Sachs Solution 

The line element of Kantowski-Saclis spacetime is 

ds^ = dt" - A\t)dr^ - B\t)\de^ - sinHd(P% (55) 

where A and B are cosmic scale factors. The corresponding Ricci scalar is 
given by 

,A 2B 2Ab 1 5^ 

'^"-^ii + ir + w + BJ + i^i. (56) 

where dot denotes derivative with respect to t. The vacuum field equations 
for Kantowski-Sachs spacetime are given by 

2B 2AB AF F ^ 
-^ + ^ + AF-F='^ ^''^ 

a b ab 1 b^ bf f _ 
-a-b^ab^b^^W^^bf'f-^- ^^^^ 

Here we also use Eq. (JT71) so that the field equations take the form 

2B 
-- + 2n 

{n + 1)— + {n^ -2n-l 
B 



B^ 


BF 

n 

BF 




(59) 


)^ 


BF 
BF 


F 1 
^F ~ B^' 


(60) 


BF 




1 


(61) 


BF 


in - 


-1)52' 



Adding these, we get 

B , ,^^2 

h (n + 1) — 

B ^ ' B^ 

Using Eq.( l23i) . we obtain one differential equation with one unknown, 
which gives 

= -^^^ + (^_i)(2n + 3)' ^^^^ 
where C3 is an integration constant. Hence the solution becomes 

ds" = {— 1 )dT^ - T^'^dR^ - T\de^ + sin2 ^02). (64) 

r4"+6 + (n-l)(2n+3) 
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where 5 = T, r = R, ^ = O and = 

The average Hubble parameter becomes here 



,n + 2 C3 1 1 

^=(^)[j^+ (^_l)(2^ + 3)r2 ]-- (65) 

while the volume scale factor turns out to be same as for the Bianchi types 
/ and /// spacetimes. The expansion scalar 9 is given by 

« = (" + 2)l^ + ^^^-^^J^^^l^ (66) 

while the shear scalar a turns out to be 

The function of Ricci scalar becomes 

while Ricci scalar takes the form 

C3(n^ + 6n + 5) 1 3n^ + 2n - 2 
^ = 2[ V^^2n^ + n-3'^- 

For a special case when n = — 5, f{R) turns out to be 

/(i?) = J_|-i?|V2. 



6 Concluding Remarks 

The main purpose of this paper is to study some Bianchi type cosmological 
models in metric f{R) gravity. Wc find exact solution of the vacuum field 
equations for Bianchi type /, III and Kantowski-Sachs spacetimes. Initially, 
the field equations look comphcated but lead to a solution using some as- 
sumptions. The first assumption is that the expansion scalar 9 is proportional 
to the shear scalar a. It gives A — S", where A, B are the metric coefficients 
and n is an arbitrary constant. Secondly, the power law relation between F 
and a is used to find the solution. Some important cosmological physical 
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quantities for the solutions such as expansion scalar ^, shear scalar o"^ and 
average Hubble parameter are evaluated. The analysis of these parameters 
in the light of solar system constraints is not done in this paper. However, 
it would be worthwhile to check the consistency of these parameters with 
Wilkinson Microwave Anisotropy Probe (WMAP) data. Also, the solutions 
can be compared with Lematre-Tolman cosmology [33j that describes the 
inhomogeneity in the universe on many scales. 

The general function of Ricci scalar, f{R), is also constructed in each 
case. For some special cases, it is found that f{R) models include some 
root powers of the Ricci scalar. It has been shown [IQ] that the model 
f{R) = R^^'^ is consistent with cosmological results. In particular, it is 
possible to obtain flat rotation curves for galaxies and consistency with solar 
system tests. However, more work has still to be done to find the viability 
of these models and consistency with solar system. 

The models of the universe in Eqs. fl29|49|IM|) are non-singular at T = 0. 
The physical parameters if, 9 and a are all infinite at this point but the 
volume scale factor V vanishes. The general function of the Ricci scalar is 
finite while the spatial part of the metric vanish at T = 0. The expansion 
stops for n = —2 in all models. The models indicate that after a large time 
the expansion will stop completely and the universe will achieve isotropy. 

2 

The isotropy condition, i.e., ^ — )■ as T — i- oo, is also satisfied in each case. 
Thus we can conclude from these observations that the models start their 
expansion from zero volume at T = and the volume increases with the 
passage of time. 
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